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Abstract 

We study the inner-structure of a charged black-hole which is formed from 
the gravitational collapse of a self- gravitating charged scalar-field. Starting 
' with a regular spacetime, we follow the evolution through the formation of an 

o ■ 

O . apparent horizon, a Cauchy horizon and a final central singularity. We find 

o 

00 

^! The no-hair theorem, introduced by Wheeler in the early 1970s, states that the 



sl null, weak, mass-inflation singularity along the Cauchy horizon, which is a 
precursor of a strong, spacelike singularity along the r = hypersurface. 



external field of a black-hole relaxes to a Kerr-Newman field characterized solely by the 
black-hole's mass, charge and angular-momentum. This simple picture describing the ex- 
terior of a black-hole is in dramatic contrast with its interior. The singularity theorems of 
Penrose and Hawking predicts the occurrence of inevitable spacetime singularities as a 
result of a gravitational collapse in which a black-hole forms. According to the weak cosmic 
censorship conjecture 0| these spacetime singularities are hidden beneath the black-hole's 
event-horizon. However, these theorems tell us nothing about the nature of these spacetime 
singularities. In particular, the final outcome of a generic gravitational collapse is still an 
open question in general relativity. 

Until recently, our physical intuition regarding the nature of these inner singularities was 
largely based on the known static or stationary black-hole solutions: Schwarzschild {space- 
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like, strong and unavoidable central singularity), Reissner-Nordstrom and Kerr {timelike, 
strong singularity). Further insight was gained from the work of Belinsky, Khalatnikov and 
Lifshitz who found a strong oscillatory spacelike singularity. However, a new and drasti- 
cally different picture of these inner black-hole singularities has emerged in the last few years, 
according to which the CH inside charged or spinning black-holes is transformed into a null, 
weak singularity P-^. The CH singularity is weak in the sense that an infalling observer 
which hits this null singularity experiences only a finite tidal deformation Nevertheless, 
curvature scalars (namely, the Newman- Penrose Weyl scalar \E'2) diverge along the CH, a 
phenomena known as mass-inflation 0. The physical mechanism which underlies this CH 
singularity is actually quite simple: small perturbations, which are remnants of the gravita- 
tional collapse outside the collapsing object are gravitationally hlueshifled as they propagate 
in the black-hole's interior parallel to the CH |[T^ (the mass-inflation scenario itself includes 
in addition an outgoing radiation flux which irradiates the CH. This outgoing flux represents 
a portion of the ingoing radiation which is scattered inside the black- hole). 

Yet, it should be emphasized that despite of the remarkable progress in our physical 
understanding of the mner-structure of black-holes, the evidence supporting the existence 
of a null, weak CH singularity is mostly based on perturbative analysis. Thus, it is of 
interest to perform a full non-linear investigation of the inner-structure of black-holes. The 
pioneering work of Gnedin and Gnedin |[ll|] was a first step in this direction. They have 
demonstrated the existence of a central spacelike singularity deep inside a charged black- 
hole coupled to a (neutral) scalar-field. Much insight was gained from the numerical work 
of Brady and Smith who studied the non-linear evolution of a (neutral) scalar-field on 
a spherical charged black-hole. These authors established the existence of a null mass- 
inflation singularity along the CH. Furthermore, they have shown that the singular CH 
contracts to meet the central r = spacelike singularity. More recently, Burko |12[ studied 
the same model problem. His work improve the numerical results given in 0, namely, he 
found a good agreement between the non-linear numerical results and the predictions of the 
perturbative approach. 
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Despite the important results achieved in these numerical investigations the mass- 
inflation scenario has newer been demonstrated explicitly in collapsing situation. It should be 
emphasized that these numerical works begin on a singular Reissner-Nordstrom spacetime 
and the black-hole formation was not calculated there. The main goal of this paper is to 
demonstrate explicitly that mass-inflation takes place during a dynamical charged gravita- 
tional collapse. 

We consider the gravitational collapse of a self-gravitating charged scalar-field 0. While 
charged collapse is not expected in a realistic gravitational collapse, it is generally accepted 
that the similarity between the inner structure of a Reissner-Nordstrom black-hole and a 
Kerr black-hole indicates that a charged collapse might be a simple (spherical) toy model 
for the more realistic generic rotating collapse. In previous papers we have investigated 



the linear [[T3| , p!^ and non-linear evolution of a charged scalar-field outside a charged 
black-hole. The results given in these papers, and in particular the existence of oscillatory 
inverse power-law charged tails along the black-hole outer horizon suggest the occurrence of 
mass-inflation along the CH of a dynamically formed charged black-hole. Thus, this model 
is suitable to establish our main goal. 

Our scheme is based on double null coordinates. This allows us to begin with regular 
initial spacetime (at approximately past null infinity), calculate the formation of the black- 
hole's event horizon, and follow the evolution inside the black-hole all the way to the central 
singularity, which is formed during the collapse. Thus, this numerical scheme makes it 
possible to test directly (and for the first time) the conjecture that the mass-inflation scenario 
is an inevitable feature of a generic gravitational collapse. 

The physical model is described by the coupled Einstein-Maxwell-charged scalar equa- 
tions. We express the metric of a spherically symmetric spacetime in the form [|16| 



ds'^ = —a{u, v^dudv + r{u, vYdVt^ , (1) 

where m is a retarded time null coordinate and v is an advanced time null coordinate. We 
fix the axis r = to be along u = v. The remaining coordinate freedom is the freedom of 



the choice of v along some future null cone, i.e. along some fixed u = const outgoing null 
ray. It should be noted that for v ^ M our null ingoing coordinate v is proportional to the 
Eddington-Finkelstein null ingoing coordinate Ve- Following Hamade and Stewart we 



formulate the problem as a system of first-oideT coupled PDEs. To do so we define auxiliary 
variables d, f, g, s, x and y: 

d=—, f = ru, g = r^, s = V47r(t), x = s^, y = Sy . (2) 
a 



Using these variables one can generalize the neutral Hamade and Stewart scheme fl^]. In 
terms of these new variables the Einstein equations expand to 

i?l = rf.-{f-£^ + ^ + lixy* + x*y) + Ueaisy* - s*y) = , (3) 
^2-r/. + /, + la^-^ = (4) 

and 

E3 = g,- 2dg + ry*y = . (5) 
The electromagnetic potential a{u, v) = Aq is given by the Maxwell equations 

Ml ^ a. - 1^ = , (6) 
where the charge q{u, v) is given by 

M2 = qy-ier^{s*y-sy*) =0 . (7) 
The Hawking mass function m{u, v) is given by 



m = -{1 + ^ + ^r^r^) . (8) 



Finally, the wave-equation for the charged scalar-field becomes 



5*1 = ryu + fy + gx + ieary + ieags + —a^qs = . (9) 
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The definition-equation (||) yields 

Dl = d- — = 0, D2 = g -r^ = 0, D3 = y - = . (10) 
a 

The initial conditions include the specification of y{0, v) and d{0, v) along an outgoing 
u = null ray. We assume d{0, v) = 0, which fixes the remaining freedom in the coordinates. 
The boundary conditions on the axis r = {u = v) are g = —f = ^a, x = y,a = q = and 
ar = Sr = (on axis). 

The evolution of the quantities d and y are determined by Eqs. El and SI, respectively. 
We then integrate Eq. Dl outwards from the axis along an outgoing [u = const) null ray 
to find a. Eqs. E3 and D2 are used to obtain g and r. Finally, we evaluate the quantities 
s,q,a,f and x by integrating Eqs. D3, M2, Ml, E2 and SI, respectively (x is evaluated 
from SI using the relation = yu)- The integration in the ^-direction is carried out using a 
fifth-order Runge-Kutta method, while the integrals in the f-direction are discretized using 
the three-point Simpson method p8[| . 

Figure |I] displays the radius r{u^ v) as a function of the ingoing null coordinate v along a 
sequence of outgoing {u = const) null rays. All the outgoing null rays originate from the non- 
singular axis r = 0, i.e. we start with a regular spacetime (this situation is in contrast with 
previous numerical works, where the evolution begins on a Reissner- Nordstrom spacetime). 
One can distinguish between three types of outgoing null rays in the rv plane: (i) The 
outer-most (small-n) rays, which escape to infinity, (ii) The intermediate outgoing null rays 
approach a fixed radius rcniu) at late-times f — > oo. This indicates the existence of a 
CH in these spacetimes. (iii) The inner-most (large-w) rays, which originate on the non- 
singular axis r = and terminate at the singular section of the r = hypersurface. These 
outgoing rays reach the r = singularity in a finite v, without intersecting the CH. This 
situation is in contrast with the Reissner-Nordstrom spacetime, in which it is well-known 
that all the outgoing null rays which originate inside the black-hole intersect the CH. One 
should also note that in contrast with the Reissner-Nordstrom spacetime, where the CH is a 
stationaryimW hypersurface, here rcniu) depends on the outgoing null coordinate u, i.e. the 
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CH contracts 0. This dramatic difference in the causal structure of the present collapsing 
spacetime compared with the Reissner-Nordstrom spacetime is attributed to the outgoing 
flux of energy-momentum carried by the charged scalar-field. 

To understand better the causal structure of our dynamical spacetime we display in Fig. 
^. the countor lines of r{u, v) in the VM-plane. The outer most contour line corresponds 
to r = 0, where its left section (a straight line u = v) is the non-singular axis, and its 
right section corresponds to the central singularity at r = 0. It should be emphasized that 
this central singularity forms during the gravitational collapse. The singularity at the r = 
hypersurface is clearly a spacelike one for r^ is negative along this section. The vanishing of 
indicates the existence of an apparent horizon (which is first formed at n ~ 1 for this specific 
solution). The CH itself is a null hypersurface which is located at f ^ oo. Its existence is 
indicated by the fact that the intermediate outgoing null rays (in the range 1 ^ n ^ 2.1 for 
this specific solution) terminate at a finite (n-dependent) radius rcniu). The singular CH 
contracts to meet the central (r = 0) spacelike singularity (along the n ~ 2.1 outgoing null 
ray). Thus, the null CH singularity is a precursor of the final spacelike singularity along the 
r = hypersurface 0. 

The behaviour of the mass function m{u, v) along the outgoing null rays is displayed 
in the top panel of Fig. ^. This figure establish explicitly the exponential divergence of 
the mass function (and curvature) along the CH in a dynamically collapsing spacetime. To 
our knowledge, this is the first explicit demonstration of the mass-inflation scenario in a 
collapsing situation starting from a regular spacetime. The mass function increases not only 
along the outgoing (n=const) null rays (as v increases) but also along ingoing (v=const) 
null rays (as u increases). 

The weakness of the null mass-inflation singularity was first predicted by Ori ||^,^. It is 
demonstrated in the bottom panel of Fig. |^. This figure displays the metric function g^y 
along an outgoing null ray, where V is a. Kruskal-like ingoing null coordinate. Clearly, guv 
approaches a finite value as the CH is approached {V 0). This confirms the analytical 
analysis of Ori, according to which a suitable coordinate transformation can produce a 
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non-singular metric. 



Figure |^ displays the Ricci curvature scalar (on the axis): 



R{u, v) 
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- [Re{y*x) + ealm{s*y)] , 



(11) 



a 



as a function of T, the proper time on axis WT^ : 



T{u) = / a{w, w)d' 







(12) 



The curvature on the axis diverges in a finite proper time. Thus, the initially regular axis is 
replaced by a (spacelike) singularity along the r = hypersurface. 

In summary, we have studied the gravitational collapse of a self-gravitating charged 
scalar-field. We calculate the formation of an apparent horizon, followed by a weak, null, 
mass-inflation singularity along the contracting CH, which precedes a strong, spacelike sin- 
gularity along the r = hypersurface. Our results give a first explicit confirmation of the 
mass-inflation scenario in a dynamical collapse that begins with regular initial conditions. 
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FIG. 1. Null rays in the rv-plane. One can distinguish between three types of outgoing null 
rays: The outer-most, which escape to infinity, the inner-most which terminate at the singular 
section of the r = hypersurface and the intermediate outgoing null rays which approach a 
(u-dependent) finite radius, indicating the existence of a CH. All the null rays originate from the 
non-singular axis r = 0. 
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FIG. 2. Contour lines of the coordinate r in the fu-plane. The r = contour hne is indicated 
by a thicker curve. Its left section (n = v) represents the non-singular axis, while its right section 
corresponds the the central spacelike singularity. The apparent horizon is indicated by the vanishing 
of Tv The (singular) CH (a null hypersurface, located at v — > oo, and indicated by the approach 
of outgoing null-rays to finite values of r) contracts to meet the central spacelike singularity (in a 
finite proper time). 
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FIG. 3. The CH singularity. The top panel displays ln(m) vs. advanced time v, along a 
sequence of outgoing null rays. The exponential growth of the mass- function demonstrates the 
appearance of the mass-inflation scenario The bottom panel displays the metric function guV 
along an outgoing null ray. The finite value approached by the metric functions is in agreement 
with the simplified model of Ori 0,^, and demonstrate the weakness of the null mass-inflation 
singularity. 
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FIG. 4. The Ricci curvature scalar R{u,u) as a function of the proper time on axis. The 
divergence of the curvature on the axis indicates that the initially regular axis is replaced in a finite 
proper time by a central (spacelike) singularity. 
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